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We develop a first-principle equation of state of salt-free polyelectrolyte solution in the 
limit of infinitely long flexible polymer chains in the framework of a field-theoretical 
formalism beyond the linear Debye-Hueckel theory and predict a liquid-liquid phase 
separation induced by a strong correlation attraction. As a reference system we 
choose a set of two subsystems - charged macromolecules immersed in a structure¬ 
less oppositely charged background created by counterions (polymer one component 
plasma) and counterions immersed in oppositely charged background created by poly¬ 
mer chains (hard-core one component plasma). We calculate the excess free energy 
of polymer one component plasma in the framework of Modified Random Phase Ap¬ 
proximation, whereas a contribution of charge densities fluctuations of neutralizing 
backgrounds we evaluate at the level of Gaussian approximation. We show that our 
theory is in a very good agreement with the results of Monte-Carlo and MD simu¬ 
lations for critical parameters of liquid-liquid phase separation and osmotic pressure 
in a wide range of monomer concentration above the critical point, respectively. 
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I. INTRODUCTION 


It is well known that thermodynamic and structural properties of polyelectrolyte in ho¬ 
mogeneous solutions are mainly determined by the correlation attraction of like-charged par¬ 
ticles which are related to long-range Coulomb interaction^I^®SIIIHni Moreover the collective 
effects in semi-dilute regime of polyelectrolyte solutions resulting from strong fluctuations of 
monomer concentration, due to the chain structure, play a crucial rol d^*^°l The necessity to 
take into account the two effects simultaneously complicates the theoretical description of 
the thermodynamics of polyelectrolyte solution in a wide range of monomer concentration 
and temperature. 

It is well known that in salt-free flexible chain polyelectrolyte solutions a liquid-liquid 
phase separation due to strong correlation attraction can take place. This effect has been 
described in several theoretical workj^ ^^ * ^^ * ^^* ^^ and has also been shown in Monte-Carlo 
simulations^. However, to the best of our knowledge theoretical models correctly predicting 
critical parameters of such a purely electrostatics driven phase transition are still missing. 
The reason why theoretical predictions of the critical parameters signihcantly deviate from 
results of Monte-Carlo simulation are three-fold. 

Firstly, in a series of studies the contribution of the correlation attraction to the total 
Helmholtz free energy has been omittecP^^^®® or described at the level of Debye-Hueckel 
approximation?2^i2l2SI2S]^ although the critical point of the electrostatic phase transition is 
located far from the range of applicability of a linear Debye-Hueckel theorj^. In worlP^ 
the contribution of the correlation attraction was described at the level of mean-spherical 
approximation (MSA) that allowed to obtain a good agreement with simulation results for 
critical temperature but underestimated the critical monomer concentration. Secondly, most 
theoretical approacheJ^^^^^^^^ffiSMl ignore the effects of the monomer concentration fluctua¬ 
tions in the regime of semi-dilute polymer solution. However, these effects, as was already 
mentioned above, should be important for solutions of polyelectrolytes as well as for neu¬ 
tral polymer^niSI jn order to evaluate the contribution of the non-electrostatic interactions 
to the total free energy of the solution the authors usually resort to some mean-held type 
approximationwhere the huctuation effects are fully ignorecP^^. Finally, all the 
above mentioned theoretical models are based on the implicit assumption that contributions 
of electrostatic and non-electrostatic interactions are independent. We would like to stress 
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that a methodology based on a ’’construction” of the total free energy using expressions ob¬ 
tained from different theoretical approaches cannot be strictly justihed from hrst principles 
of statistical thermodynamics. 

In referencd^ addressed an interesting observation that one can obtain reasonable values 
of the critical parameters if one uses the modihed random phase approximation (MRPA) for 
the calculation of the electrostatic contribution into the free energy of the solution. MRPA 
was developed in work^^SED contrast to the standard random phase approximation 

(RPA) contains a concept of ultraviolet cut-off in the procedure of integration over vectors 
of reciprocal space. The parameter of ultraviolet cut-off determined from the condition of 
equality between number of degrees of freedom of the system and number of collective vari¬ 
ables which contribute into the total free energy enumerated by elements of reciprocal space. 
Therefore, the cut-off parameter within MRPA is related to the number density of particles. 
It is surprising that MRPA allows one to obtain very accurate interpolation formulas for 
excess free energy of one component plasma (OCP) and hard-core one component plasma 
(HCOCP) in a wide range of temperatures and number densities - from regime of linear 
Debye-Hueckel theory to the strong coupling regimd^^*^. 

In workP^ in the framework of a field-theoretical approach the interpolation formulas for 
excess free energy, correlation length, and effective Kuhn length of the segment of polymer 
chain for semi-dilute solution of neutral polymer chains were obtained. It should be empha¬ 
sized, that the interpolation formulas obtained by authors contain all limiting laws predicted 
in earlier works within the scaling approach®^. Using the widely used Carnahan-Starling in¬ 
terpolation formula for the free energy of a hard sphere system and the interpolation formula 
for excess free energy of semi-dilute polymer solution mentioned abovd^ offer now new op¬ 
portunities for theoretical description of more complicated polymer systems in the framework 
of thermodynamic perturbation theory (TPT). 

The presented observation in referencd^ and results proposed in referenc^^ motivated 
our attempt to develop a statistical theory of salt-free polyelectrolyte solutions which would 
simultaneously take into account the effect of correlation attractions of charged monomers 
and counterions beyond the linear Debye-Hueckel theory as well as the collective effects 
related to the monomer concentration fluctuations in the semi-dilute regime. 

We formulate the theoretical model based on a variant of TPT formalism, which has been 
developed in work^^^*^. As the reference system we have chosen a set of two independent 
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subsystems - charged polymer chains immersed in a structureless oppositely charged back¬ 
ground created by counterions (polymer one component plasma) and counterions (which 
we model as charged hard spheres) immersed in oppositely charged background created by 
polymer chains (hard-core one component plasma). We calculate the excess free energy of 
the polymer one component plasma in the framework of MRPA. We show that our theory is 
in a very good agreement with Monte-Carlo simulations results reproducing critical parame¬ 
ters of the electrostatic phase transition and MD simulation results reproducing the osmotic 
pressure as a function of monomer concentration in a region that is above the critical point, 
respectively. 

II. THEORY 

A. Hamiltonian 

In this section we present the theoretical model based on the TPT formalism for the 
thermodynamics of a salt-free flexible chain polyelectrolyte solution. The system contains 
charged polymer chain and counterions, where we assume that the charged polymer chain is 
immersed in a good solvent modelled by a structureless dielectric medium with the dielectric 
permittivity £. 

Moreover, we take into account an excluded volume of monomers and counterions mod¬ 
eling them as hard spheres of diameter d. For simplicity we also assume that each monomer 
has a fixed charge —e, whereas the counterions carry an opposite charge -t-e ( monovalent 
counterions). Thus the free energy of the solution within the canonical ensemble is 

F = -kBT\nZ, ( 1 ) 

where the partition function Z takes the form 

with 

H = + //W + i + i (ft Vft) + (ftKft) 

being the Hamiltonian of the system, where the following short-hand notations have been 
introduced 

{fVcg) = j dx j dy/(x)W(x-y)^(y), (3) 
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K(x-y) = 


£|x-y| 


is the Coulomb potential, 


N ^ 

Pp(x) = -e^ / cis(5(x - Ri(s)) 
i=i i 


(4) 


(5) 


is the microscopic charge density of the monomers, 


Nc 

Pc(x) = e^(5(x- 
i=i 


( 6 ) 


is the microscopic charge density of the counterions. 


E / / *1*2^ (R.(i>i) - R,(i>2)) 


Nr, 


NN 




(7) 


0 0 


is the Hamiltonian of the monomer-monomer excluded volume interaction; w = 27id^/3 is 
the second virial coefficient of the excluded volume interaction of hard spheres of diameter 
d] P = l/ksT is a reciprocal temperature, ks is Boltzmann constant. 




( 8 ) 


is the Hamiltonian of counterion-counterion excluded volume interaction; 

oo, |r| < d 


f4,(r) = 


0, |r| > d. 


(9) 


is the hard-core potential. 


N 




( 10 ) 


0 


is the Hamiltonian of counterion-monomer excluded volume interaction. Moreover, integra¬ 
tion measure over the conhgurations of polymer chains has been introduced: 


Np N 

drp(-)= VR,.. VRjv^e ^=^0 ( 


( 11 ) 


where b is the Kuhn length of the segment, Np is the number of polymer chains with N as 
the degree of polymerization; 

j dr,(.) = a;’"' / *^ / *».(■) (12) 

V V 
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is the integration measure over the phase space of counterions, Nc is the total number of 


counterions, in the system of volume V. The symbol J 'DR(..) in (11) denotes a functional 
integration over the polymer chain configurations. Moreover, the following normalization 
condition must hold 


VUe 


N 

^ f dsR'^(s) 


= V. 


(13) 


In order to take into account a monomer-counterion excluded volume interaction we use an 
assumption 


dr, / dTpe-^^-(■) ^ / < / dr;(-). 


where 


and 


Vm.r. — Vr. 


<(■) = A 


-3Afc 


^ drl.. j 

Vf Vf 


dT'J-) = / T>Ri.. / VRme 


N Np 


(14) 


(15) 


(16) 


= Vc = Vm = 'n'd^/6. The approximation (14) is a simplest way to include into 
the theory depletion effect!^. The symbol J' in (16) denotes that an integration performed 
over the free volume for counterions Vf = V — NcVmc- In the renormalized measure J drp(-) 
we have the following normalization condition 


N 

— f dsR^(s) 

VRe ” « =Vf. 


(17) 


B. Construction of Perturbation Theory of System 

Adopting the notation POCP for the polymer one component plasma and HCOCP for 
the hard-core one component plasma, the model Hamiltonian (|^ can be rewritten as 

H = Hpocp + HhCOCP + Hpert, (18) 

where 

Hpocp = + - (ppVcPp) + [ppVcpc) + - iPcVcPc) (19) 

is a Hamiltonian of polymer one component plasma (POCP), 

HHCOCP = + - (pcVcPc) + (pcycpp) + 2 (pp^cPp) ( 20 ) 
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is a Hamiltonian of hard-core one component plasma (HCOCP); 


Hpert PcS^c^ Pp) 

is a perturbation part of the total Hamiltonian; pc = eNc/Vf is an average charge density of 
counterions, Vf = V — N^Vcm is a free volume for counterions (monomers), pp = —pc is an 
average charge density of monomers, 6pp(x.) = Ppipi.) — pp is a local charge density fluctuation 
of monomers, Spd'x.) = Pc(x) — pc is a local charge density fluctuation of counterions. We 
would like to stress that a set of charged flexible polymer chains immersed in a structure¬ 
less neutralizing background we shall call throughout the paper as polymer one component 
plasma (POCP). Using the formalism of TPT we shall use a set of two independent ideal 
subsystems as a reference system - POCP of charged macromolecules and HCOCP of coun¬ 
terions. Using the TPT formalism for the two independent ideal subsystems as a reference 
system, we obtain 

Z = (22) 

where symbol ((■))/? denotes an averaging over microstates of the reference systemP^l 

{(•))« = (23) 

The partition function of the reference system is then given by 

Zr = exp [—(3Fpocp — I^Fhcocp] ■ (24) 

Hr = HHCOCP + Hpocp is a Hamiltonian of the reference system. Thus an expression for 
the total Helmholtz free energy F within the TPT approach has a following form 

^F = f 3 Fpocp + f^Fncocp - In ^ . (25) 


C. Free Energy of Reference System 


The expression for density of excess free energy of HCOCP has a following forirP^ 


^HCOCP, ex 


At] — 3 

(1 — pY 4 


In (0 cZqPc) - cPc 



3 / cPcZq 

2 V 0 


3/2 

arctan 


0 


cTr.Zn 



(26) 
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where Tc = Ib is a plasma parameter, rj = Tid^n'^/Q is packing fraction of coun¬ 
terions, n'^ = Urn/ g”"” ^ is a renormalized concentration of monomers (counterions), 

Ib = e^/ekBT is a Bjerrum length, Um = NpN/V is concentration of monomers (counteri¬ 
ons), and c = 2/3 (4/7r^)^'^^ is a number constant. Moreover, the following notations have 
been introduced 

(1-# 


Zn = 


0 = 1 + 


1 + 4?7 + 4:7]^ — + 7]'^ ’ 

61b ~ 7])*{16 — Ht] + 

5d (1 + 4?7 + 47/2 — 47]^ + 7]‘^y 


(27) 


(28) 


The equation (26) very accurate describes the thermodynamic properties of HCOCP in wide 


ranges of number density and temperature - from range of applicability of Debye-Hueckel 
theory to strong coupling regimd^. Applying an analogous approach to calculation of free 
energy of POCP as well as for HCOCI®, we arrive at the following expression (see Appendix 
A): 


POCP,ex — d-POCP,e.v. + d-POCP,cor, 


(29) 


where the correlation contribution into the free energy of POCP can be determined by the 
expression 

^POCPyCor ’^mk'Bd'^ CT^ 1 ( 30 ) 


where 


= I (^ - In (l + a; + ^) ) + 
3 a /2 00 (\/ooa — a / cjct ^ — 4 ) 




ooa^ + \/ua^ (ooa^ — 4 ) 

arctan t /--h 

2oo 


3a/2 00 (^/uja + Voj(T^^-4) Iuia"^ — J (cjcr^ — 4) 

: arctan *' 


a^/^-y/cjcr — ^Joo (cncr^ — 4/ 


2oo 


(31) 


1/3 


00 = cPm, Pm = Tc = /b ^ j ig a plasma parameter of monomers, o = (97r^n(„)^/^y^47r62Q;/9 
is a dimensionless parameter. The non-electrostatic part d^pocp,e.v. of fhe excess free energy 
of POCP can be determined by the expression which obtained by Edwards and Muthukumar 
within a field-theoretical approach® 

J 2 


-T, 


POCP,e.v. 

knT 


2471^3 9)77 


9wn'^ ^ wn'^ 


( 32 ) 


























where correlation length ^ takes the form 
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2 \ 47rnL6^Q! 


+ 


9 


AnriLb'^a 


+ 


2ab‘^ 

SwnL 


( 33 ) 


and expansion factor a of neutral polymer chain satisfies the following equation 


= 




(34) 


The system of coupled equations (32), (33) and (34) determines the density of free energy 


of solution of flexible polymer chains with excluded volume and qualitatively describes ther¬ 
modynamic properties of the solution in the case of good solvent from regime of semi-dilute 
polymer solution to concentrated regim^^. 


D. Fluctuation Correction to Free Energy of Reference System 

Further, let us calculate the fluctuation correction to the free energy of the reference sys¬ 
tem at the level of the Gaussian approximation. Using the standard Hubbard-Stratonovich 
transformation we rewrite the partition function of the solution in the following form 

■PT 


Z = Z 


R 


c 




where 


U(x-y) = 


0 14(x-y) 

K (x - y) 0 


is a vector of auxiliary fields. 


T(x) = 


5p(x) = 


TUx) 


6pp{x) 

(5pc(x) 


(35) 


(36) 


(37) 


(38) 


is a vector of local charge density fluctuations. Moreover, the additional notations have been 
introduced 

r f „ .. .X 


and 


a, 5 

(T,(5p) = j cix^T„(x)5p„(x). 


(40) 
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An operator V ^ can be determined by the following integral relation 


[ ^ “ y) = '^« 7 ^ (x - y) • 


( 41 ) 


Performing some cumbersome calculations (see Appendix B) we arrive at the expression 


ex — POCP,ex + HCOCP, ex “ 1 “ peril 


(42) 


where the perturbation part of the total density of free energy of the solution at the level of 
Gaussian approximation takes the form 


f^^pert n 


dk 

W)- 


In 1 - ^ ((5p,(k)(5p,(-k))popp {(5ft(k)4pc(-k)) 


HCOCP 


K^(k) , (43) 


where 14(k) = ^ is a Fourier image of Coulomb potential, (5/5p(k)5/5p(—k)) 


POCP 


IS a 


Fourier image of correlation function of microscopic charge density fluctuations of POCP, 
and {6pc{k)Spc{—k))fjcocp ^ Fourier image of correlation function of microscopic charge 
density fluctuations of HCOCP. The latter has been obtained within the Gaussian approxi¬ 
mation in workpS ^nd can be written in the following form 

Zoe^n'^k^ 


(5p.(k)5pe(-k)) 

1 


(44) 


^HCOCP- j^Zo' 

The corresponding expression for POCP also obtained within Gaussian approximation (see 


Appendix A) has a following form 

{6p,{k)6p,{-k)) 


POCP 


n'^e^Sik)k^ 

k? + Aim'llsSik) ’ 


(45) 


where 5'(k) = is a well known expression for the structure factor of semi-dilute 

polymer solution?^. 


Substituting (44) and (45) into (43) and calculating integral we arrive at 

1 


pert 


12<3 


(G(m, v) — G{u, 0)), 


(46) 


where u = \/AttI v = AtiIbu'^Zq^, and function G(m, v) has a form 


G(m, v) = 1 - 1 -n^-l 


■\/2 {2v? -|- (n^ -f 1) (n^ — 1)^ — Au^ — — 1 j sj+ 1 + — ly — 4^ 


Ayju^ + 

\p2 — (n^ -I- 1)\J (n^ — 1)^ — Au^ — — 1 j + 1 — — ly — 4^ 


( 47 ) 
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E. Equation of State 


The expression (42) in combination with expressions (26 29) and (46 47) dehnes the excess 


free energy of salt-free polyelectrolyte solution in the limit of infinitely long polymer chains 
{N —)■ cx)) which taking into account the correlation attraction of charged particles, monomer 
concentration fluctuations in the regime of semi-dilute solution and mutual influence of the 
latter effects. It should be noted that the contributions of electrostatic and non-electrostatic 
interactions are not independent. The osmotic pressure can be obtained by standard way 


n — ^id + ^exi 


(48) 


where 

Hid = n'^ksT 

is an ideal part of the osmotic pressure and 


H-ex — TT'm 




dn. 


-Te- 


(49) 


(50) 


m / j' 


is an excess osmotic pressure of the solution. 


III. NUMERICAL RESULTS AND DISCUSSION 

Turning to the numerical calculation we introduce the dimensionless monomer concentra¬ 
tion hm = temperature T = ksTed/e'^, and the second virial coefficient w = wd~^. In 

addition we assume that b = d. Moreover, we determine the dimensionless osmotic pressure 

n = nedVeT 

Fig.l shows the osmotic pressure as a function of monomer concentration for two tempera¬ 
tures T. At sufficiently low temperatures the appearance of a Van-der-Waals loop indicates a 
liquid-liquid phase separation in the solution. Using the equality between chemical potentials 
and the osmotic pressures of coexisting phases we obtain a coexistence curve of the liquid- 
liquid phase separation with an upper critical point (see Fig.2). The obtained coexistence 
curve exhibits a highly asymmetric shape, which is in agreement with Monte-Carlo computer 
simulation results b}^. Fig.2 also shows for comparison the values of the critical points ob¬ 
tained within different theories (including the present theory) and Monte-Carlo computer 
simulations. The values of the critical parameters obtained by extrapolation of simulation 
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results for —)■ oo in the workP^ are: = 0.2 and hm,c = 0.09. The Muthukumar theorj^ 
sufficiently underestimates both critical parameters: = l/Odvr 0.005. The the¬ 

ory developed by Jiang and co-author^ based on MSA, as has already been mentioned in 
the Introduction predicts the critical temperature in close agreement with simulation results 
but underestimates the critical monomer concentration: Tc = 0.245 and nm,c = 0.01. The 
present theory yields the critical parameters Tc = 0.2 and hm,c = 0.12 which are in close 
agreement with Monte-Carlo simulation results. 

It is instructive to regard how our theory can describe the thermodynamic properties of 
polyelectrolyte solution in the region above the critical point. Fig. 3 shows the comparison 
between the dependencies of the osmotic pressure on monomer concentration calculated 
within present theory and within MD simulation at T = 1.2 in wide range of monomer 
concentratiorP^. As seen from Fig. 3, our theory gives very good agreement with the results 
of MD simulation. Also our results are in very good agreement with a numerical calculation 
within PRISM theorjP^. It should be noted, that two qualitatively different regimes of the 
osmotic pressure behavior take place. At a region of the small monomer concentration where 
the entropy of the mobile counterions gives a basic contribution into the osmotic pressure 
there is a linear dependence on monomer concentration (IT ~ u-m), whereas in concentrated 
regime we observe an essentially nonlinear behavior (H ~ hm^)- The predicted scaling law 
has been also obtained within a Monte-Carlo simulatioiP^ and numerical calculations based 
on PRISM theorjP^. 

IV. SUMMARY 

We have presented a first-principle statistical theory of salt-free flexible chain polyelec¬ 
trolyte solutions in the regime of a good solvent. The approach is based on the thermody¬ 
namic perturbation theory formalism and predicts a liquid-liquid phase separation arising 
from strong correlation attraction of charged particles. Taking into account the correlation 
attraction of charged particles beyond the linear Debye-Hueckel theory and strongly cor¬ 
related monomer concentration fluctuations in the regime of semi-dilute polymer solution 
we have developed an equation of state for the salt-free polyelectrolyte solution that accu¬ 
rately predicts the critical parameters of the liquid-liquid phase separation and the osmotic 
pressure in a wide range of monomer concentrations in the region above the critical point. 
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We would like to stress that our theory provides accurate values of the critical parameters 
without the concept of charge renormalization (e.g. counterion condensation). As is well 
known, the counterion condensation is an essentially nonlinear effectP^® that cannot be 
described within the linear Debye-Hueckel type of theory. However, our theoretical model 
going beyond the simple linear theory, implicitly takes into account the effect of charge 
renormalization. 

A limitation of the presented theoretical model is related to the fact that it can be applied 
only to polyelectrolyte solution of sufficiently long polymer chains (i.e., in a case N —)■ oo). 
However, as was shown in theoretical workd^^®®® and conhrmed by Monte-Carlo and MD 
simulations in reference^^QEiEl degree of polymerization of the polyelectrolyte chains 
has only a weak effect on the thermodynamic properties of salt-free solution. In other 
words, if the degree of polymerization is higher than ~ 30, the thermodynamic quantities 
of polyelectrolyte solution in the regime of semi-dilute solution only weakly depend on the 
length of the polymer chains. 

The obtained accurate equation of state can be applied to compute thermodynamic prop¬ 
erties of polyelectrolyte solutions with a low concentration of low-molecular weight salts 
using the ionic strength as a small parameter of perturbation within the thermodynamic 
perturbation theorj^. In addition, the obtained equation of state can be used to describe 
inhomogeneous polyelectrolyte solutions within the classical density functional theory. Re¬ 
placing the equation of state of the simple one-component plasma!^ by the equation of 
state into the free energy functional instead of allowing to adequately take into account the 
contribution of the correlation attraction. 

In conclusion we would like to speculate about the reasons why a set of two independent 
one-component plasmas can be used for theoretical description of thermodynamics of salt- 
free polyelectrolyte solution near the critical point. As was predicted in the referencd^ and 
conhrmed by Monte-Carlo simulation in workP^, due to the strong correlation attraction the 
charged macromolecules near the critical point form a network of polyelectrolyte aggregates. 
Thus, in this case the charged network can be regarded as a weakly huctuating background 
neutralizing the charge of counterions, whereas the mobile counterions immersed in the 
polymer network can be regarded as a neutralizing background for charged macromolecules. 
Thus, in order to evaluate the contribution of the remaining effects that are related to the 
huctuations of charge densities of the neutralizing backgrounds which undergo only slight 
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fluctuations the Gaussian approximation can be applied. 

In conclusion we would like to estimate the critical parameters in physical units for the 
real experimental systems. We take the following parametes for the polyelectrolyte and 
counterions 6 ~ d = 0.25 nm which approximately corresponding to sodium polystyrene- 
sulfonate. As an example of polar solvents, where the discussed liquid-liquid phase transition 
can take place, we take N,N-dimethylformamide {e ~ 36.7) and acetonitrile {e ~ 37.5). For 
N,N-dimethylformamide we obtain the following value of the critical temperature ~ 
364 K. For acetonitrile we get ~ 356 K. The critical concentration in both cases 
fim,c — 13 M. Thus this purely electrostatic liquid-liquid phase transition can be realized at 
experimentally accessible conditions. 
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V. APPENDIX A: DERIVATION OF FORMULAS (30) AND (45) 


In this Appendix we derive the expressions for the excess free energy (30) and the corre¬ 


lation function of the charge density fluctuations (45) of POCP. 
Starting from the partition function of POCP 


^POCP 


— f>~f^^POCP,e.v.—0HpoCP,el 


(51) 


where 


H 


ksTw 


Nr, 


N N 


POCP,e.v. 


^ dsidsrS (Ri{si) - Rjis;)) 


(52) 


*0-10 0 

is the monomer-monomer excluded volume interaction part of the total Hamiltonian, 

Np N 

f f f E fdsR^(s) 

J dPp(-) = J PRi..y VR^pC -1" (•) 

is the probability measure over configurations of the Gaussian polymer chains, 

HpQCP^ei 2 ^Pp^cPp) “1“ “1“ 2 (P^^bcPfo) 


(53) 


(54) 
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is the electrostatic interaction part of the total Hamiltonian due to monomer-monomer, 
monomer-background and background-background interaction, pi, = — (/5p(x))^^^ is the 
charge density of the neutralizing background, f3Eseif = ^ 2 I^t ^(Q) electrostatic 

self-energy of monomers and 


K(x - y) = 


- y| 


(55) 


is the Coulomb potential. 


We now turn to the calculate the partition function (51) in the framework of TPT for¬ 


malism which was proposed in the referenc^^ choosing as a reference system a set of neutral 
polymer chains with excluded volume interactions. As announced in the main text of the 
article we are making use of the interpolation formula for the free energy of such system 
that takes into account the effect of strongly correlated monomer concentration fluctuations 
in regime of semi-dilute polymer solution!^. 

Hence, we obtain 


^POCP 


^n.ch.f' 




n.ch. 


(56) 


where 5pp{x.) = Ppipi.) — (pp(x))^^^ is the local charge density fluctuation of the monomers; 
Zn.ch. is the partition function of a solution of neutral polymer chains. The symbol ((■))nc/i 
denotes averaging over microstates of solution of neutral polymer chains. 


Applying to (56) the standard Hubbard-Stratonovich transformation we arrive at 
Zpocp = I , 

where 


(57) 


C= / IS (>”''• ■'") (58) 

is a normalization constant. For sake of simplicity the following short-hand notations 


and 


V) = <iyv>(x)V ‘(x-y)»’(y) 


{Sppip) ^ / dxSpp{x)ip{x). 


(59) 


(60) 


have been introduced. The kernel of the reciprocal operator 1/ ^ can be determined from 


the relation 


cizW(x-z)V; (z - y) = 5(x - y). 


(61) 
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Truncating the cumulant expansion in (57) at the second order 


=exp 


^ ~ 2^"^ (^Pp^pp)c^) - 


we obtain 


ZpoCP ~ Zn.chfi 


/3Es^if [ ^p)-^{p{SppSpp)^ip) 


C 


( 62 ) 


(63) 


where the symbol {(■))p denotes the cumulant averagd^. Performing the calculation of the 


Gaussian integral (63), we arrive at 
Zpocp = Zn.ch. exp 


V 

J 


dk 


(27r)3 


k: 


k: 


d-iMi + d(i + W) 


k2 


(64) 


where Kp = ^'''^Tdpn^. We shall use the following approximation for the structure factor of 
the semi-dilute polymer solution 

9 1 


k2 47ra62k2’ 

where the following approximate expression for correlation radius has been used 


(65) 


AnrimOib'^ 


( 66 ) 


The evaluation (66) follows from the equation (33) for the case of the semi-dilute polymer 
solution. The expansion factor a of the neutral polymer chain in the case of semi-dilute 
polymer solution as follows from (34) can be evaluated as a ~ (dS/vr)^'^^ , so that 

evaluation of the correlation radius takes the form 

1/4 




/ 243 \ 
Vl0247r2 ) 


{wb 




(67) 


that is in agreement with the well-known scaling result!®. 

Further, following the idea proposed in the reference^®® introduce an ultraviolet cut¬ 


off in the integral over the wave vectors k in (64). In order to obtain the cut-off parameter 
A, we assume that total number of collective variables dpp{k) which contribute to the total 
free energy of the POCP is equal to the total number of degrees of freedom of POCP, i.e. 
3NNp. Therefore, we arrive at the following relation 

A 


2V 


dvr 


dkP = SNN., 


( 68 ) 
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from which we obtain the expression for cut-off parameter 


A = ( 97 rSi„)‘''’. 


( 69 ) 


The prefactor 2 on the left-hand side of (68) originates from the summation (integration) 
of both the real and the imaginary parts of the collective variable 6pp{k.). Thus, using the 


ultraviolet cut-off (68) and taking the integral in (64) we obtain 


^POCP — ^POCP,e.v. + ^POCP,cor) 

where 

^POCP,cor (J (cj, (t) 

is a correlation contribution and 

9i^, = I (^ - In (l + a; + ^) ) + 

3V2 OJ (a/w(T - a/cuct^ - 4) 


(70) 


(71) 


cr^/^-^CJcr -|- y/uT (cncr^ — 4^ 


arctan 


-|- y/cjcr^ (cna^ — 4) 

2uj 


3a/2 oj iy/ua -I- \/ooa^ — 4) ju:a^ — 

: arctan \ /-- 

2u 


^ _ .. (72) 

(T^/2 uja — \Juj (cncr^ — 4) 

uj = cTm with Tra = Ip {^^TTUm being the plasma parameter of the monomers, c = 
2/3 (4/7r^)^'^^ a number constant and a = (97r^nm)^^^\/47ra6^/9 is a dimensionless parameter. 


Turning to the calculation of the correlation function (45) of the local charge density of 


POCP at the level of the Gaussian approximation we employ the method of the generating 
functionaP^. The partition function of POCP immersed in an imaginary auxiliary external 
held iV’(x) takes the form 


^Eseif / „-§iSppVcSpp)+i{-tp5pp) 


Z[^] = Zn.ch.e^ 

After introducing the following generating functional 


n.ch. 


Z[0]’ 


(73) 


(74) 


with Z[0] = Zpocp the correlation function of local charge density of POCP can be deter¬ 
mined from the double functional derivative of the generating functional ^['0], he. 


{6pp{x.)6pp{y)) 


POCP 


5'0(x)(5'0(y) 


P=o 


(75) 
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Using the Hubbard-Stratonovich transform of the partition function Z[tlj\ we obtain 




( 76 ) 


Thus, we obtain the following functional representation of the generating functional 

^ f ^^i{^p5pp)+ii^5pp)'^ 

Zpert J C 


(77) 


where 

'^pert j ^ 

is the perturbation part of the partition function. Performing a shift of the integration 




variable ip ^ cp — ip we obtain 

f ^g-^uv-V)+^(</^i4-V)-2)3uvt-V) ^giM/3p)^ 

Zpert J C 

Using the cumulant expansion in the integrand 

^gi((/7(5pp)^ _ gi/(ix¥>(x){<5pp(x))^-| /dx/dy{5pp(x)5pp(y)}^(p(x)v5(y)+.. 

and truncating at the second order we obtain 


(79) 


(80) 


(;[p] ^ ^ f V)+^(vV7 V) 


Zpert .1 C 


(81) 


The kernel of the integral operator D ^ can be determined in the following way: 


D (x-y) = U^ (x-y)+/3((5pp(x)%(y))^, 


(82) 


Analogously, truncating the cumulant expansion at the second order in Zpert and taking 
the Gaussian integral we obtain 


Q[ijj] ~ exp 


■^(vv-'v) + ^(vv 


-^DV~^ 

c ^ ^ c 


Ip) 


= exp 


-^(V-GVO 


(83) 


where 


G = K"' - V-^DV-^ = V-^ (/ - DV-^) = V-^ (j - {V-^ + A)-^ UU') = 

= UU' (/-(/ + VeA)-^) = UU' (/ + KA)-' (/ + KA -I) = 

= V-^ (J + KA)-' KA = VaXA (/ + KA)-' = A (/ + KA)-', (84) 
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and A is an integral operator which can be represented as 


A(x -y)= (3 {5pp{yi)5pp{y))^ = (3e^nmS{yi - y), (85) 

I is the identity operator, S'(x — y) is the structure factor of a neutral polymer chain in 

i(k) 


good solvent. Using the relation (84) we obtain 




,ik(x-y) 


( 86 ) 


where 

i(k) = /3e2n^5(k). (87) 

Thus we arrive at the Fourier transformed correlation function of the charge density 
fluctuations of POCP 

e^nm5'(k) 


(%(k)5pp(-k)) 


POCP 


l + ^e%„V;(k)S(k)' 


( 88 ) 


Substituting the Uc(k) into (88) we arrive at (45). 


VI. APPENDIX B: DERIVATION OF FORMULA (43) 


In this Appendix we present the derivation of the formula (43) for the perturbation part 


of the total free energy. The partition function of the solution within TPT has the form 


Z = Z 


R 


c 




(89) 


Further we have introduced the following definition for the perturbative part of the partition 
function 


Qpert ry 




c 


(90) 


where the identity (5p(x))^ = 0, and the definition of the inverse operator W 


-1 


(X - y) = U (X - y) + /5 (Spa(x)Sp^(y)} 


R 


(91) 


have been taken into account. Further, calculating the Gaussian integral (|90|) we obtain 

V r dk ^ det(U-Hk)) 


Qpert exp 


In 


(27r)3 det(lU-i(k)) 


(92) 
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with the matrices W ^(k) and V ^(k) being 


w-^ (k) = 


/3(%(k)%(-k)) 
1 

Vc(k) 


v-^ (k) = 


POCP Vc(k) 

/?(5pe(k)<5p,(-k)) 

HCOCP, 


0 


V-c(k) 


.Vt(k) 


0 


(93) 


(94) 


Further, calculating the determinants of the matrices (93) and (94) and substituting results 
into (92) with (3J^pert = —^InQpert, we obtain the formula (43). 
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13049. 



Figure 1. The osmotic pressure IT as a function of the monomer concentration at different tem¬ 
peratures T. At sufficiently low temperatures the presence of a Van-der-Waals loop indicates a 
liquid-liquid phase separation. 
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Figure 2. The coexistence curves of salt-free polyelectrolyte solution and critical points that are 
obtained within present theory and Monte-Carlo simulations^ and critical point calculated within 
theory of Jiang et.aipZl. The critical point calculated within the present theory is in very good 
agreement with Monte-Carlo simulations results. 



Figure 3. Osmotic pressure 11 versus monomer concentration hm as obtained by the present theory 
MD computer simulation^, and the PRISM theorj^ at T = 1.2. 
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